The Nernst effect from fluctuating pairs in the pseudogap phase 
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The observation of a large Nernst signal in cuprates above the superconducting transition tem- 
perature has attracted much attention. A potential explanation is that it originates from super- 
conducting fluctuations. Although the Nernst signal is indeed consistent with gaussian fluctuations 
for overdoped cuprates, gaussian theory fails to describe the temperature dependence seen for un- 
derdoped cuprates. Here, we consider the vertex correction to gaussian theory resulting from the 
pseudogap. This yields a Nernst signal in good agreement with the data. 
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For most of the doping phase diagram, high temper- 
ature superconductivity in the cuprates emerges from a 
normal state where an energy gap is already present. 1 
The origin of this 'pseudogap' is the subject of much 
debate. 2 One of the most interesting observations in 
the pseudogap phase is the existence of a large Nernst 
signal. 3 ' 4 The Nernst effect is the generation of a trans- 
verse electric field by a thermal gradient in the presence 
of a magnetic field perpendicular to both. Since vor- 
tices carry entropy, it is natural to attribute such a large 
Nernst signal in proximity to the superconducting tran- 
sition temperature T c to vortex-like excitations. 3 ' 4 More- 
over, invoking vortices is consistent with the Nernst sig- 
nal smoothly going through T c . On the other hand, it is 
not clear whether vortices give an adequate description 
of the physics of fluctuating superconductors except very 
near T c . 5 The free vortex density above the Kosterlitz- 
Thouless temperature increases exponentially with tem- 
perature, a dependence which is inconsistent with the 
near power-law decrease of the actual Nernst signal above 
T c . Moreover, the recent observation of a negative Nernst 
signal for underdoped YBa2Cu306+ y has further compli- 
cated the story. 6 This negative signal has been argued to 
be a consequence of density wave reconstruction of the 
Fermi surface. 

Here, we take the point of view that the dominant con- 
tribution to the Nernst signal in the pseudogap phase is 
indeed due to fluctuating pairs. This is supported by the 
close correspondence of the Nernst signal with the fluctu- 
ational diagmagnetism. 7 On the other hand, we note that 
although existing theories based on Ginzburg-Landau or 
diagrammatic approaches 8-10 give a good description of 
the Nernst data for overdoped cuprates, they do less well 
for underdoped cuprates. We attribute this to the pres- 
ence of the pseudogap. 

As discussed in Ref. 10, it is the direct contribu- 
tion from fluctuating pairs - the Aslamazov-Larkin (AL) 
contribution 11 - which governs the Nernst signal over a 
wide range of temperatures above T c , so we focus on that. 
The AL contribution to the Nernst coefficient is obtained 
from the electric current-heat current Kubo response ker- 
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The latter can be expressed in terms of cor- 



graphs) connected by interaction lines (pair fluctuations) . 
The vertex block can be expressed as Tr^^'^B} with 
B = (vGGG) where v is the Fermi velocity, and (...) indi- 
cates disorder averaging. The factor 7 differentiates the 
electric vertex — e) from the heat vertex 7^, which 
we will discuss below. Disorder averaging in B leads to 
the presence of two Cooperons and the renormalization 
of the free electron Greens function G by impurity scat- 
tering. In order to account for the pseudogap, we replace 
this G which was previously used to compute this block 
by the broadened BCS Greens function 



G(k,e) 
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as this gives a good description of photoemission data in 
the pseudogap phase. 12-14 In Eq. (1), A*, is the momen- 
tum dependent pseudogap and e = e + T sgn(e) with T 
the scattering rate. By recomputing the electromagnetic 
vertex block with this G, we find that B is renormalized 
by a function of A/T where A is the maximum value 
of the pseudogap. Assuming a T independent A and 
T ~ T as observed in photoemission, 12,15 this renormal- 
ization results in a fluctuation Nernst signal which drops 
off considerably faster with temperature than the gaus- 
sian result. As we show, this gives a good description of 
the Nernst data for underdoped cuprates. 

We assume the standard expression for the pair prop- 
agator whose retarded component has the form 16 
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responding electric and heat vertex blocks (triangular 



Here, e = ln(T/T c ), N is the density of states, and 
7] = irD/8T where D is the diffusion constant. The 
Nernst coefficient can be expressed in terms of the electri- 
cal (tr) and thermoelectric (d) tensors as v = (a xy a xx — 
a xx a xy )/H(al x + cr 2 y ) ~ a xy /Ha xx . The second (ap- 
proximate) expression, which becomes exact if particle- 
hole symmetry is present, gives a good approximation 
to v for the case of superconducting fluctuations since 
a xy 3> a xx and a xx 3> o xy (see Ref. 16 for a corre- 
sponding discussion). The transverse thermoelectric co- 
efficient, a xy — a xy + cM z /T, consists of two indepen- 
dent contributions: the response of the total current to 
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the applied electric and magnetic fields (a xy ), and the 
magnetization currents as derived from the equilibrium 
magnetization M z . Like Refs. 9,10, we focus on the first 
contribution 



H 
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assuming here weak field limit. The electric current-heat 
current Kubo response kernel 

A xy (Q, iQ m ) = -Ae 2 T ^ B x (q)B 2 {q){ioj n + ifl m /2) 
[L(q - Q x ,iw n )L(q, iu> n )L(q, 

+ L(q, iuj n )L(q, iu> n + ifl m )L(q + Q x ,iuj n + iO TO )](4) 

is written in Matsubara representation, where we have 
assumed that the heat vertex is —iuo n /{2e) times the 
electric vertex. 17 In Eq. (4), the Greens function block 
B, whose rcnormalization is the subject of this paper, 
is assumed to be independent of frequency. This ap- 
proximation is formally exact in the immediate vicin- 
ity of the transition temperature. Nevertheless, a rigor- 
ous extension of this approximation to a wider range of 
temperatures above T c demonstrates that the Ginzburg- 
Landau result, v Ah <~ (T — T c )~ , remains valid even 
far from T c if one substitutes T — T c by the more gen- 
eral expression Tln(T/T c ). We note that in the gaussian 
approximation 
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For the following linear response calculation we expand 
the propagators in Eq. (4) to the leading order in Q, 
namely L(q ± Q x ,iuj n ) — > ±Qd q!c L(q,iuj n ), and noting 
that 



dL(q,iw n ) 
dq x 



= -B^(q)L 2 (q,iLU r , 
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we thus find from Eq. (4) 



A xy (Q,in m )= (7) 
-4e 2 QT B x (q)B^(q)B 2 y (q)(ioj n + iQ m /2) 

[L 3 (q, iu> n )L(q, iuj n + iQ m ) - L(q, iuj n )L 3 (q, iu> n + iCl m )] 

Performing the summation over the Matsubara fre- 
quency uj n by using contour integration, T J2u n f(^ w n) = 
^iq § c duj coth ^/(w) with two branch-cuts at Im(aj) = 
and Im(u;) = — f2 m , followed by an analytic contin- 
uation il m — > and keeping only the linear in 51 
contribution from A xy (Q, f2), one finds for the transverse 
thermoelectric coefficient 
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After the remaining momentum and energy integrations, 
we get the result (restoring K) 
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where £h = \Jhc/eH is the magnetic length and £q L = 
r\l In (T/T c ). For B = B G , this is the well known expres- 
sion for the Nernst effect from fluctuating pairs. 9 ' 10 

At this point, all we have done is to rederive the gaus- 
sian expression. The reason we have done this is to 
demonstrate explicitly where the Bi current vertices en- 
ter. We now discuss the renormalization of Bi due to the 
pseudogap. The expression for the vertex is 11 

Bi(q, to, fi) = T ^ «i(p)G(p, e)G(p, s + Sl)G(q -p,u-e) 

p,e 

(10) 

where s is the fermionic loop frequency, u> the bosonic 
frequency that enters the fluctuation propagator, D, the 
external field frequency (set to zero for the dc response), 
and G the Greens function. 18 We will now recalculate 
Bi using the pseudogap Greens function. Substituting 
Eq. (1) in the previous equation, taking the dc limit, and 
approximating G(q — p,uj — s) w G(q — p, —e) we obtain 
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Keeping the term to linear order in q, 
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Converting the p sum to an integral, we have 
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where A# = Acos2i9 (i.e., we assume a d-wave pseu- 
dogap). Next we perform the £ integral by introducing 
£ = fix with fi — \Js 2 + A'^ 



Bi(q) ~ ^N Tv 2 Fqi J2 J 
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The x integral is trivial, and we find 
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and F/A = y/3T/T* with the coefficient chosen so that the 20 40 6 ° so 100 o 20 40^60 so 100 

spectral gap disappears at T* . 



The angular integral is easily performed, leading to 
Bi(q) ■ 

[(4e 2 + 2A 2 )E(X) - e 2 K(X)] (16) 
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where A = A/\/e 2 + A 2 and K and E are elliptic func- 
tions. As the sum converges at n > A/T, one can ap- 
proximate the elliptic functions by their value at zero 
argument, which is ir/2. We will also take the 'zero T' 
limit by converting the Matsubara sum to an integral, 
noting that the dominant T dependence comes from the 
T dependence of T. This gives 
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12 "Jo e 2 (e 2 +A 2 ) 3 / 2 
and after the remaining integration results in 
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where we have exploited the BCS relation £0 = Hvp/nA. 
Inserting the gaussian expression for B, we obtain 
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In the limit of small r <C A, the term in parenthesis 
reduces to 2A/T. Since V ~ T and 77 ~ 1/T, then in 
this limit, the ratio is a constant, and one obtains the 
same functional form for the Nernst as in the gaussian 
approximation. On the other hand, as T increases, the 
ratio decreases from unity as can be seen in Fig. 1. This 
leads to a Nernst signal which decays more rapidly than 
the gaussian result, since three renormalized B vertices 
enter the expression for the Nernst: 
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FIG. 2: Nernst signal, v, versus T for L^-zSr^CuCU for four 
different values of the doping, x. 4 The curve for x=0.17 is a fit 
to the gaussian expression for a xy . The other curves include 
the vertex correction as described in this paper. 



We now consider the Nernst data for La2- £C Sr :E Cu04. 4 
The advantage of these data is that after the normal car- 
rier background has been subtracted, the Nernst signal 
is positive, and therefore complications due to density 
wave reconstruction can to first order be ignored. In 
Fig. 2, we show the Nernst signal, v = —E y /(H Z V X T), 
for four different dopings. For the overdoped compound, 
the gaussian expression for a xy fits the data quite well, 
but not for the three underdoped compounds where the 
pseudogap is present. Instead, we find that the corrected 
expression provides a good description of the data (with 
pseudogap temperatures T* ranging from about 200 to 
300 K). 

We now turn to a brief discussion of the bosonic con- 
tribution to the conductivity. As with the Nernst, the 
paraconductivity observed in underdoped cuprates well 
above T c falls off like T~ s with 5 « 3. 19 We recall that 
when extended to the higher temperature regime, gaus- 
sian theory (in 2D) predicts for the Aslamazov-Larkin 20 
contribution to the conductivity a xx = e 2 /8n 2 ln 3 (T/T c ) 
while for the Maki- Thompson 21 contribution, o\^, T = 
7r 2 e 2 ln(l/7 ¥ ,)/1921n 2 (T/T c ) where 7^ is the dephasing 
rate. In either case, the ln _Tl (T/T c ) decay is too slow to 
explain the experimental data. We now argue that the 
same vertex renormalization can account for the faster 
decay of the fluctuational conductivity. We start from 
the definition of the conductivity in the linear response 
regime a xx = — Im[A^.(n)]/0 where the current-current 
response kernel is 

A xx (itt m ) = -4e 2 T ^ b I(q) l (Q^ iu n )L(q, iu„ + ifl m ) 

(21) 
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After Matsubara summation and analytic continuation, 
this reduces to 

o p 2 r+oo 

= —E B 'W J ducoth—lmL R (q,uj) 



: [L R (q, u) + ST) + L A {q,uj - Si)] 
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In the dc limit, we expand L R (q,u! + Sl) + LA.(q,ui — Sl) ~ 
2iSld u lmL R (q, u), integrate by parts over w, and as a 
result find for the AL contribution to the conductivity 

p 1 r+oo j 

°**=^ B2 ^L^i ilmL « M]2 (23) 

This is formally the same expression as in gaussian theory 
except for the vertex function now defined by Eq. (18). 
We thus find as a result (restoring h) 
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where the renormalization factor (B / B G ) 2 provides a 
faster power-law decay, consistent with the data. 19 

We would like to conclude with the observation that al- 
though we assumed a 'BCS' expression for the pseudogap 
Greens function, G, any theory of the pseudogap with 
a T independent d-wave like gap and a scattering rate 
proportional to T will yield results equivalent to those 
derived here. Similar conclusions have been reached in 
regards to the fermionic contribution to various transport 
properties in the pseudogap phase. 22 

In summary, we note that the Nernst signal and fluc- 
tuational conductivity for underdoped compounds drops 
more rapidly with temperature than predicted from a 
gaussian theory of fluctuating pairs. This discrepancy 
is nicely accounted for by a vertex correction to the 
fermionic current block due to the pseudogap. 
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